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^rj ' It has been proposed that the adiabatic loading of a Bose-Einstein Condensate (BEC) into an opti- 

cal lattice via the Mott-insulator transition can be used to initialize a quantum computer [D. Jaksch, 
et al, Phys. Rev. Lett. 81, 3108 (1998)]. The loading of a BEC into the lattice without causing band 
excitation is readily achievable; however, unless one switches on an optical lattice very slowly, the 
optical lattice causes a phase to accumulate across the condensate. We show analytically and nu- 
merically that a cancellation of this effect is possible by adjusting the harmonic trap force- const ant 
O |' of the magnetic trap appropriately, thereby facilitating quick loading of an optical lattice for quan- 

tum computing purposes. A simple analytical theory is developed for a non-stationary BEC in a 
harmonic trap. 



I. INTRODUCTION 



o 

Experimental advances in manipulating and controlling Bose-Einstein Condensates (BECs) of dilute atomic gases 
has resulted in a remarkable series of experiments [Q . One theoretical proposal for quantum computing using atoms 
as qubits is to first load the atoms that are in a BEC into an optical lattice. Then, by varying the intensity of a laser 
used to form an optical lattice the BEC will undergo a quantum phase transition from its BEC-like superfluid state 
to a Mott-insulator state |2) . This has recently led to a seminal experiment by Bloch and collaborators [|| . 

In principle, starting with a BEC in a trap and turning on an optical lattice of sufficient well depth in a sufficiently 
adiabatic manner will prepare the Mott-insulator state. In practice, it is easy to turn on the optical lattice adiabatically 
with respect to band excitation (excitation from one band to another); however, it is substantially more difficult to 
turn on the optical lattice adiabatically with respect to quasi-momentum excitation. The second, more stringent 
form of adiabaticity requires that the optical lattice be switched on slowly with respect to mean-field interactions and 
tunneling dynamics between optical lattice sites, and hence typically requires milliseconds Q]. We will refer to the 
first form of adiabaticity as 'mterband adiabaticity' and the second form as 'mfraband adiabaticity'. The intraband 
adiabaticity condition has been demonstrated in one-dimensional lattices by Orzel et al. Q and ultimately led to the 
pioneering experimental demonstration of the Mott-insulator transition || . When not otherwise specified, the terms 
■ adiabatic and nonadiabatic in this article will refer to intraband adiabaticity. 

The goal of the present paper is to present a simple strategy for remaining in the adiabatic regime while switching 
on the optical lattice much faster than the millisecond time scales ordinarily required for intraband adiabaticity. The 
strategy is to counterbalance the switching on of the optical lattice with an appropriate change in the force constant 
O |. of the trap. This strategy is shown to correct and prevent much of the quasi-momentum excitation and resulting 
phase damage that arises from the nonadiabatic nature of the switching. 

More specifically, the switching on of an optical lattice potential can divide a BEC into many individual pieces where 
phase coherence is maintained across the whole condensate. This phase coherence can be seen by instantaneously 
dropping the lattice and looking at the momentum distribution through time of flight measurements. However, 
because of a spatially dependent change in the density and thus the mean-field per well site, one can end up with a 
quadratic phase dependence developing along the lattice direction if one does not load the lattice adiabatically with 
respect to quasi-momentum excitations Q. Elsewhere Q it has been shown, using optimal control methods, that 
one can control the phase evolution to obtain a flat phase at some final time by time varying the harmonic trap 
force-constant of a confining external (typically magnetic) trap. Here we show analytically and numerically that a 
complete cancellation of the phase development is possible by appropriately adjusting the external trap. 

This paper will focus solely on one-dimensional lattices, considering only the dynamics of the BEC along the lattice, 
and will ignore effects transverse to the lattice. It should be noted that the effects of transverse excitation will show-up 
on time scales inversely proportional to loj_ the transverse trapping frequency which is typically long compared to 
the times in the present paper. Work is now in progress toward further extending these results to two- and three- 
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dimensions. It is expected that the squeezing of the BEC into the transverse directions can also be treated using 
the above method, namely by an appropriate adjustment of the trap in those directions. 

There have been a number of recent publications of both experimental || 0, [l0| and theoretical ]?], studies 
involving the loading of BECs in one-dimensional lattices, and the resulting dynamics. This paper is related to these 
publications but focuses explicitly on a means of quickly loading an optical lattice from a BEC for quantum computing 
purposes, as well as for improving experimental signal to noise in short time experimental studies of BECs. Note that 
we consider the regime where the density of the condensate is sufficiently large that mean-field effects are not enti rely 
negligible. Experiments can be carried out in the truly dilute gas regime where mean-field effects are negligible [ p.0| . 
However, reducing the condensate density to such low values would have to be carried out adiabatically, adversely 
affecting the time to load the optical lattice from the initial (dense) BEC. 

The outline of the paper is as follows: in section [n] we define the problem. In Sec. Ill A, a simple analytical theory is 



developed for a nonstationary ID BEC in a harmonic trap. It is shown that a change in the density of the condensate 
induces a time- varying phase across the condensate that can be eliminated by a change in the harmonic force constant 



of the trap. In section HI B it is shown that the effect of switching on the optical lattice is to generate a new effective 
normalization of the BEC and an analytical expression is obtained for the modified harmonic trap force-constant 
that compensates for the new effective normalization. The analytical theory is in excellent agreement with numerical 
simulations. A modified version of the theory in the regime where the n online ar interaction is strong and hence the 
the width of the condensate differs from well to well is developed in Sec. Ill Cj . Section IV contains the conclusion. 



II. DESCRIPTION OF PROBLEM 

We consider a ID BEC confined by a harmonic trap and governed by the Gross-Pitaevskii equation 

iftj^) = (K + V + NUo\iP\ 2 M, (1) 

ft 2 r) 2 

where K — —-^-§^2 is the kinetic energy operator, V is the external potential energy operator to be discussed shortly 
and NUq is the nonlinear atom-atom interaction strength, A~ being the number of atoms and Uq = ATraoh 2 /m is the 
atom-atom interaction strength that is proportional to the s-wave scattering length a^. The BEC is initially in the 
ground state of the trap potential and is therefore stationary. An optical lattice is then switched on, having the effect 
of separating the BEC wave packet into a series of localized pieces. The potential energy operator therefore takes 
the form V(x,t) = ^mut 2 x 2 + S(t)Vo cos 2 (kx), where u>t is the trap frequency - which may be time dependent, k 
is the laser field wave number, Vq is the lattice intensity and S(t) is the function that switches-on the laser for the 
optical lattice and goes from S = at the beginning of the ramp-on of the optical potential to S = 1 at the end 
of the switching on time St s - In applications to quantum computing, one often wants to create an optical lattice 
with one atom per lattice site which will serve as quantum bits. However, due to the nonlinearity of the equations, 
the condensate wave function develops a phase that varies from lattice site to lattice site when the optical lattice 
is not turned on adiabatically Such a wave function can be represented by a superposition of quasi- momentum 
states, and a superposition of quasi-momentum corresponds to a higher energy state and thus cannot give rise to the 
Mott-insulator state. The problem we address is the elimination of this phase profile by adjusting the trap strength. 
In the coming section we analyze the evolution of BEC wave functions in harmonic traps, and consider the effect of 
switching on the optical lattice. Finally, a closed form for the precise time dependence of the trap strength that will 
insure a flat phase for the wavefunction for all times after the optical potential is fully turned on is derived. 

First, however, we transform the NLSE to dimensionless units t — > t/to, x — > x/xq and ip — ► y/xQtp where for 
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convenience we choose to = -W 1 ■ Performing these transformations we end up with a dimensionless NLSE 

i^(x,t)= (~\^ + K(t)x 2 + S(t)V cos 2 (kx) + U\^\ 2 ^, (2) 

where the trap force-constant K = w 2 tg, the field intensity V = Voto/U and the nonlinear coefficient U = NUoto/xoh, 
such that all space, time and energy quantities are now expressed in units of x$, t and h/t respectively. 1 



1 We do not determine, at this point, any specific choice of xq. Note however, that choosing xo = X/tt, the optical wave length, yields 

h 2 k 2 

2m 



for the energy units, h/tg = n k = E r which is just the recoil energy. 
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III. ANALYTICAL THEORY 



A. Dynamics of a Thomas-Fermi BEC in an Harmonic trap 

Consider a normalized Thomas-Fermi type BEC wave function in a harmonic potential of the form 



4{x,t) = \ V&V 1 -^-* ("M^ 1 , (3) 
{ (x/w) 2 > 1 

where the width w(t) and phase components b(t) and c(t) are all assumed to be time dependent. We wish to analytically 
describe the evolution of this wave function in a harmonic trap with trap force-constant K (we first consider the case 
where K is constant in time, but the equations of motion for w(t), b(t) and c(t) remain valid even if K varies with 
time). Inserting the above wave function into the dimensionless NLSE, we obtain, by considering separately the real 
and imaginary parts, two equations involving the three parameters w(t), b(t) and c(t). The imaginary part yields 

w , 2x 2 s b , 2x 2 s 
-5^(1 -^ = -0(1 -7^) = W 



where a = y 1 — , and from the real part we get 

'X 



bx 2 -C = -^r+tV+^ + Pf (1-^) 

Aw a Aw 



-^(l + =^) + (b 2 + K)x 2 + U-^(l-^). (5) 

4ur w Aw ur 

In going to the last line we expanded \ in a Taylor series in a;/u>, truncating after the second order. Comparing 
separately the coefficients of x° and x 2 , we obtain the following two equations of motion for b(t) and c(t): 



I 3U 

Aw 2 Aw 



(7) 



Taking a time derivative of Eq. (Q) and using Eq. (^|) we find 



bw + bw 
_L _3C/ 

2 W 3 4y;2 

(8) 



Kw 



dw 

with the effective potential V e (w) defined as 

V e (w) = --^ + ^ + ]-Kw 2 . (9) 
4w 2 Aw 2 

The time evolution of the wave function width, w, can therefore be easily determined by considering the form of the 
potential V e (w). Furthermore, by defining 

p = wb, (10) 

we can formulate the equations for the conjugate variables w and p as a Hamiltonian system of equations with 
H(w,p) = p 2 /2 + V e {w) such that 

w = ^H=p (11) 

d d 
^ dw dw e 
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FIG. 1: a) Ve(ui) with stationary point wo and examples of turning points wi and W2 marked. The asymptotic curves correspond 
to the contributions of the two dominant terms in V e and highlight the way changes in these terms effect the dynamics, b) V e ' 
chosen such that w' = wi by adjusting K. c) V e ' obtained from (a) as a consequence of change in wave function normalization, 
d) a new V e ' obtained from (c) by also changing K to compensate for the change affected by the normalization change shown 
in (c). 

Consider now the potential V e (w) in (||) plotted as curve (a) in Fig. [j]. The potential consists of a well centered 
around the stable point wq ~ (Jtf) 5 * This can be most easily obtained by setting 6 = and b = in Eq. (||) and 
solving for w while noticing that the first term on the RHS of (|^) is small compared to the rest and can therefore 
be neglected. With initial wave function width u>(0) = wq, where wq is the width of the Thomas- Fermi stationary 
ground state of the trap, the wave function will remain stationary throughout. However, if the initial width equals 
some other value, an oscillatory motion of w round the stationary point wq will develop. The phase curvature b will 
also oscillate with w obtaining its maximum value when w(t) — wq and vanishing when w approaches its turning 
points, wi and W2- 

If an abrupt change in the trap force-constant can be made, K — > K' , at the exact point in time when b(t) = 0, 
i.e., when w(t) is at one of its turning points, e.g., w — wi, then it is possible to change the potential so as to freeze 
the flat phased wave function and make it stationary. This can be obtained by choosing K' — such that wi is 

the stationary point of the new potential V' R (curve (b) in Fig. |l]). 

Another scenario to be considered is the following. We begin with a stationary flat phased wave function residing 
at the stationary point wq of the potential. Imagine now the hypothetical possibility of abruptly changing the 
normalization of the BEC wave function from unity to n. This would be equivalent to a change in the potential 
V e — > V e ' affected by changing U — > nil. It is obvious that this change will shift the stationary point to some new 
value w' = 3 (see curve (c) in Fig. [I]) and that the wave function currently positioned at wq will no longer be 

stationary under the new potential. In order to compensate for this change and keep the wave function stationary 
one can adjust the trap force-constant and set K' = nK such that ratio U/K remains constant and the stationary 
point w' = wq will not shift (see curve (d) in Fig. |l|). 

We show in the following section that turning on an optical lattice corresponds to a change in the normalization 
of the wave function, so that the above scenario corresponds precisely to our goal of achieving a flat phased BEC 
loading of an optical lattice. It should be noted that the above analysis ignores gravity which can be assumed to be 
orthogonal to the lattice direction. However, even if gravity is along the lattice direction a similar analysis holds but 
requires an additional linear offset. 
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B. Switching on the Optical Lattice 



Quickly switching on the optical lattice causes the BEC wave function, which initially has a Thomas- Fermi form of 
an inverted harmonic potential, to split into a series of localized pieces each residing in a lattice well. As the overall 
normalization of the wave function must remain unity, the displaced population from areas between the lattice wells 
builds up within the wells such that the density in these regions increases dramatically (see Fig. ||). However, if we 
neglect the local lattice structure and consider solely the global nature of the BEC wave function, we see that it 
retains its quadratic shape, and the change in the wave function brought about by the existence of the optical lattice 
can be viewed as a stretching of the Thomas-Fermi wave function in the vertical direction (see Fig. ||) . This picture 
is based on a separation of scales in the spatial dimension which is a consequence of the fact that the l eng th of each 
lattice well, ^ = r , is much smaller than the scale of the total wave packet, w (see for example Ref. JllJ). It is for 
this reason that we can treat first the local structure of the wave function in each well and then consider separately 
the overall global evolution of the wave function. 

The idea is therefore to view the wave function on a level coarser than the lattice site dimension, averaging out 
the local lattice structure of the wave function. This procedure yields a new Thomas-Fermi type wave function 
I V'giob | 2 = (|V'| 2 )ioc = fS(l — iu 7 ) differing from the initial one by a modified normalization factor n (see Fig. ||). The 
evolution of this wave function can then be analyzed using the results of the previous section. 

This procedure can also be viewed as a spatial-averaging out of the local structure of the Hamiltonian operator 



(#>ioc = (T + lattice + V t + U\iP\ 2 ) loc = T av + V^ tice + V t + [/<|Vf>i oc . The harmonic trap potential V t is constant 



on the local scale and is therefore unaffected by the averaging. If the average kinetic and lattice potential energies 
per particle, T av and V^tticc: are constant from well to well, these contributions to the energy can be absorbed into 
the chemical potential /i, resulting in just the averaged global mean-field playing-off, on the global scale, against the 
trap potential as in a simple Thomas-Fermi procedure. The trap must then be adjusted to compensate only for the 
varying mean-field across the BEC wave function. 

In obtaining this simplified picture we distinguish between two opposite scenarios occurring on the local scale. In 
many cases, when considering the dynamics along the direction of the one-dimensional lattice, the mean-field within 
each well is negligible in comparison with the kinetic and potential energies along this direction. This occurs for tight 
optical wells, e.g., short wavelength and strong intensity such that \/Vk 2 U/w, where w is the width of the BEC. 
The local wave function can then be well approximated by a Gaussian with a "well-independent" width implying 
that the locally averaged kinetic and lattice potential energies are also "well-independent" . In carrying out the above 
procedure we then find that the global wave function is a stretched image of the initial one, as described above. 

In the opposite regime the mean-field within each well can no longer be neglected. In these cases the calculations 
are more involved and do not yield the simplified picture presented here of a mere stretching of the wave function. 
Instead, a distortion occurs which must be treated explicitly. We therefore delay discussion of this scenario and 
provide a more general treatment in the next section. 

In the following we wish to determine the normalization factor n in terms of the optical lattice parameters V and 
k. Consider the initial Thomas-Fermi wave function \i(j\ 2 = xr(l — The number of atoms in the region of each 
lattice well determined by its position Xi is 



Assuming that the local population becomes trapped in the well during the switching on of the optical lattice, we 
can then consider the local normalization factor per well as constant throughout the evolution. Assuming too that 
the wave function at each lattice site is localized after the optical lattice has been switched on, we can ascribe to each 
lattice site a local wave function, (f>i(x), which is normalized to fy. In order to obtain an average norm per well we 
define the local probability function Pi{x) which is just the local wave function normalized to unity 



Averaging out the local structure using the local probability function, Pi(x), we obtain the coarse-grained wave 
function, 



Jj(Xi) = \^{xi)\ 2 - 



Awk 




(13) 



Pi(x) = -|<Mx)| 2 . 



(14) 




(15) 
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FIG. 2: Analysis of the BEC wave function in an optical lattice, ipi n it and ipfinai are the wave functions before and after 
applying the optical lattice, and 4>i(x) is the local wave function within a specific well. (Gaussian approximation) %j) s \ b is the 
global Thomas-Fermi type wave function after averaging out the local details. 



Note, the limits of integration in the above integral should be restricted to a single well but due to the gaussian-like 
nature of the wavefunction <f>i(x) the specific limits are unimportant. Note that in evaluating the integral, r)i was 
considered constant as it is only slowly varying on the local scale. 

In many cases the local wave function can be well approximated by a Gaussian 



<l>i(*) = ] /j^e-^B L j*, (16) 

where A is the width and the wave function normalizes to the local normalization factor r\i (see inset in Fig. |J). A 
is of typically on the order of but smaller than A <C w and is therefore small compared with the width of the total 
wave function, so T]{Xi) is only slowly varying with respect to x and can be considered constant within any given 
lattice site. Averaging out the local structure we obtain the coarse wave function ipgiob(%i) which we now show to be 
of Thomas-Fermi type 



|Vgiob(^)l 2 = (l^l 2 ) l0 c = J K{x)M*)\*d* 
rn 1 f _ 2 (^f 



/2ttA 

3 /F 1 x? 



3n x 2 
Aw w 



( 17 ) 



In going from the third to the fourth line we used the explicit form of r\i given in (|i3|). Comparing the last two lines 
we find the modified normalization to be 



7T 1 

2 fcA' 



(18) 



It remains to determine the local width A of the wave function within each lattice site in terms of the external 
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parameters. It can be shown analytically (see appendix |X|) that 



where W(x) is the Lambert W function p2|, so that the normalization factor n is finally given by 



n=U ~\ . (20) 

The effect of switching on the optical lattice on the dynamics of the wave function can now be viewed as changing 
the normalization of the initial wave function from unity to n. 





FIG. 3: Evolution of the wave function (amplitude and phase) as a consequence of switching on the optical lattice. Note the 
development of a quadratic phase profile. 

If the switching-on time St s is short compared to the global nonlinear time scale t^ l |L3) so as no substantial phase 
evolution occurs during this time, the transformation of the normalization constant can be considered abrupt and the 
dynamics of the wave function parameters w(t) and b(t) are raised from the initial potential curve V e to V e ' (curves (a) 
and (c) respectively in Fig. |lj) by the change U — > nU as described in the previous section. If no further adjustments 
are made, the wave function will begin to evolve on the potential curve V' e and develop phase as seen in Fig. ^[ In 
order to cancel this effect one can compensate for the change of normalization by adjusting the trap force-constant 
to K — nKo (curve (d) of Fig. Q). In Fig. ^ we plot the switching-on function of the optical lattice and the change of 
the trap force-constant AK = K — Kq as a function of time. The evolution of the wave packet under this sequence 
of events is plotted in Fig. |5|, from which it is evident that the phase remains constant throughout the evolution for 
the correct tuning of the trap force-constant. 

In the simulations presented here we have taken N — 1.5 x 10 6 sodium atoms, a scattering length of ao = 2.8 
nm, and a trap of average frequency 59.26 Hz. Using these values the Thomas-Fermi approximation to the chemical 
potential, /xtf, can be calculated and the nonlinear interaction time becomes t^L = Ji/fJ-TF = 96.2 /is. In order to 
preserve the time scales in the ID model as they are in 3D reality, we follow Ref. and replace the nonlinear 
coefficient NUq by C^tf^tf, where the Thomas-Fermi radius xtf = \/ 2/xtf / mutf gives the size of the condensate 
and the factor C carries the dependence of the simulation on the dimensions and is for our ID case C — r ( 2 +f/2) = f 
& 

We take the optical lattice wavelength to be A = 589 nm and choose xq = A/7r = 2/k such that to = h/E r , 
where E r = is recoil energy. The optical lattice is switched on in a time 8t s ~ 20/zs to the final intensity of 
Vq = 10.94 E r . In units of xq = A/tt, to and h/to = E r , for space, time and energy quantities respectively, we therefore 
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FIG. 5: Evolution of the wave function (amplitude and phase) as a consequence of switching on the optical lattice and adjusting 
the compensating trap force-constant. The stationary flat phase is strikingly apparent. 



get the following unitless values; k = 2 for the optical wave number, K — {^ L ) 2 — 5.615 x 10 6 for the initial trap 
force-constant, V — 10.94 for the final field intensity and U = | j-^p 1 = 9.55 for the nonlinear interaction strength. 

Inserting these values into Eq. ( ]2C| ) yields the normalization factor n = 2.0866 such that the trap force-constant 
which we analytically predict to yield an optimally flat phase is K — 1.172 x 10~ 5 (85.5 Hz). This value is off by 
merely 2% from the empirically found optimal value of K = 1.151 x 10 -5 (84.8 Hz) which generates the evolution 
plotted in Fig. ||. Some small residual spatially varying phase structure remains. This structure is due to incomplete 
interband adiabaticity and can be reduced by increasing the switching-on time 8t s . 
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C. Nonlinear Regime 



We now return to the more complicated scenario where the local wave function has spatially varying contributions 
from the mean-field term. Various complications arise in this regime which must be solved individually. The main 
complication is due to the fact that when the mean-field is locally important it affects the width and shape of the 
local wave functions such that they differ from well to well as shown in appendix This implies that the average 
kinetic and lattice potential energies also vary from well to well, affecting the phase accumulation. 

Assuming that the local wave function can still be approximated by a Gaussian along the lattice direction (as is 
the case unless the local mean-field is larger than the kinetic energy) we c an use the results of appendix |a| to obtain 



the well- dependent width A (2^). This can be inserted back into Eq. (A2) to obtain the total local energy as a sum 
of its contributions: the kinetic energy T, the lattice potential energy .©lattice, the trap potential energy Et and the 
mean- field energy E m f. The chemical potential associated with a specific lattice site is [[l4| 

(T + ©lattice + E t + 2E m f) 



n 



(l-e- K/Xi ) + Vt(xi) + -^p-. (21) 



SAf 2 v ' " v " V2^A 



In order to keep the phase evolution constant from well to well one must adjust Vt(xi) such that it cancels all other x% 
dependencies (originating in A(xi) and rj(xi)) and thus makes fi independent of xi. Another complication arises from 
the fact that the optical lattice must be switched on adiabatically with respect to interband excitations (as stressed 
in the introduction), e.g., the switching-on time 5t s in our dimensionless units must be longer than £jy - This means 
that for experiments in which the lattice wavelength is large, the lower bound on the switching-on time becomes 
comparable to t^L and considerable phase evolution will occur during this time. 

To avoid the phase winding during the switching-on time one must make the trap frequency change gradually so 
as to compensate for the changing shape of the wave function at intermediate times. We assume as a zeroth order 
approximation that a transition of the magnetic trap from its initial to its final form using the same switching-on 
function as the optical lattice, will momentarily compensate for the changing shape of the wave function. The relevant 
parts of the potential terms in the Hamiltonian will take the following form (1- S(t))V t init ■ + S (t)V^ nal + S (t)V cos 2 (fcx). 

In Fig. ||we show a case where the mean- field is important. In this simulation we chose parameters as above, except 
the optical lattice wavelength and strength were changed to be A = 8 x 589 nm, 2 and Vq = 45.4 E r respectively so that 
the mean-field within each well is no longer negligible. With these parameters and following the above procedure we 
found the optimal trap shape to be of the form V t flnal = ©/{l-(l-(^) 2 )°- 8 } with K = 11.88, where w = x TF /x = 13.5 
is the width of the BEC in the units of xq introduced above. We turned on the new trap shape gradually, as described 
above, with a switching-on time of 5t s w 1 ms, the resulting constant flat phase can be clearly seen in Fig. ||. 
Some small residual spatially varying phase structure due to incomplete interband adiabaticity remains here too, and 
increasing the switching-on time St s will reduce the residual phase structure. 



IV. CONCLUSIONS 



The switching on of an optical lattice potential can divide a BEC into many individual pieces where phase coherence 
is maintained across the whole condensate. However, because of a spatially dependent change in the density and thus 
the mean-field per well site, one can end up with a quadratic phase dependence developing along the lattice direction 
if one does not load the lattice adiabatically. We have shown analytically and numerically that a cancellation of this 
effect is possible by appropriately adjusting the external trap. A simple analytical theory has been developed for a non- 
stationary ID BEC in a harmonic trap. It was shown that the effect of switching on the optical lattice is to generate a 
new effective normalization of the BEC, and hence a nonstationary condensate. Finally, an analytical expression was 
obtained for the modified harmonic trap force-constant that compensates for the new effective normalization. The 
analytical theory is in excellent agreement with numerical simulations. 

In real experiments more care is needed to account for the effects of evolution in the transverse directions. Work 
is now in progress toward extending these results to two and three dimensions. It is expected that the expansion 
of the BEC into the transverse directions can also be treated using the above method namely by an appropriate 



2 This can be accomplished by changing the configuration of the lasers that make the optical lattice from counter-propagating to inter- 
secting at an angle 9 such that sin(0/2) = 1/8 (i.e., 9 = 14.36 degrees). 
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FIG. 6: Evolution of the wave function (amplitude and phase) as a consequence of adiabatically switching on the optical lattice 
and adjusting the compensating trap (strength and shape). 

adjustment of the trap in those directions. We have detailed elsewhere how our quasi-lD calculations of the type we 
presented here model 3D aspects of the dynamics in cylindrically symmetric potentials 0, but this method can not 
describe radial excitations of the BEC that might arise due to the optical potential via the mean-field interaction. 
To the extent that radial excitations are not important, our method should be an adequate approximation to the 3D 
dynamics. 

It is not known how a small residual spatially varying phase will affect the Mott-insulator transition. The residual 
phase can be thought of as a phonon-like excitation that should be mapped onto the final Mott-Insulator state. 
Characterizing the nature of excitations in an inhomogeneous Mott-Insulator has not been done; however, the small 
residual excitations seen here are not expected to have a strong effect since the total energy of the system is only 
slightly above that of the ideal case. A more exact answer to this question can not be provided within the context of 
a mean-field approach and requires analysis using many-body approaches to the Mott-insulator transition. Moreover, 
no theoretical model exists that is completely appropriate in both the superfluid and Mott-Insulator regimes. 

APPENDIX A: CALCULATION OF LOCAL WAVE FUNCTION WIDTH A 

As in the text we approximate the local wave function by a Gaussian 

of width A, normalized to the local normalization factor rji and centered around Xi. We wish here, using the variational 
method, to determine the width of the Gaussian in terms of the optical lattice strength V and wave number k. 
We first compute the energy associated with 4>i as a function of A: 

E(A) = J 4>*{x) [~\^ + V t (xi) + Vsm 2 (k(x - a*)) + ^U\4>i(x)\ 2 ^ <f>i{x)dx 

n(r \ r (x-x t 2 / 1 ft 2 1 \ r x -m -) 2 

= (^g£) J e~^~ ^---^ + V t (x l ) + V S m 2 (k(x-x t )) + -U\<f )l \ 2 )e-^^dx 

Note that the trap potential, denoted Vt and the number of atoms in the region of the lattice well at position Xi, 
rj{xi), were extracted from the averaging integral since they are assumed constant on the local scale. According to 
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the variational principal, the determine the ground state function, 0j, which minimizes the energy E(A) with respect 
to A: 

An explicit solution of this equation for A is not possible in general; we therefore distinguish between several cases 
and make some simplifying assumptions. If, as is the case for short optical wavelength, the mean-field term becomes 
negligible with respect to the other energy terms, it can be neglected to obtain the following equation 

At '-" A ' = «!p- < A4 » 

The solution to this secular equation can be written in terms of the Lambert W function, y = W(x), which is defined 
as the inverse of x = ye y [Il2| , 



It can be seen in the inset of Fig 
important point to note is that in t 



A " f i w( -i7w^ (A5 » 

that this value for the width of the local wave function gives good results. An 
;his regime A is independent of the well position, implying that the lattice potential 
and kinetic energies per particle too are "well-independent" . This crucial point justifies our treatment of the global 
wave function i/'giob as a Thomas-Fermi approximation. 

For high density BECs and longer optical wavelengths the mean-field cannot be neglected and Eq. ( |A3|) must be 
numerically solved for A. It must be noted, however, that the resulting form for A(xi) will in general be "well- 
dependent" , implying that the kinetic and lattice potential energies per particle will also be "well-dependent" and 
thus contribute to the phase curvature accumulation. This must be taken into account when adjusting th e trap to 



counter the phase accumulation, within the non-negligible mean-field regime, as will be discussed in section III C . 
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